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Abstract. Within the class of stochastic cellular automata models of traffic flows, we look at the velocity 
dependent randomization variant (VDR-TCA) whose parameters take on a specific set of extreme values. 
These initial conditions lead us to the discovery of the emergence of four distinct phases. Studying the 
transitions between these phases, allows us to establish a rigorous classification based on their tempo- 
spatial behavioral characteristics. As a result from the system's complex dynamics, its flow-density relation 
exhibits a non-concave region in which forward propagating density waves are encountered. All four phases 
furthermore share the common property that moving vehicles can never increase their speed once the 
system has settled into an equilibrium. 
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1 Introduction 

In the field of traffic flow modeling, microscopic traffic 
simulation has always been regarded as a time consum- 
ing, complex process involving detailed models that de- 
scribe the behavior of individual vehicles. Approximately 
a decade ago, however, new microscopic models were being 
developed, based on the cellular automata programming 
paradigm from statistical physics. 

The seminal work done by Nagel and Schreckenberg 
in the construction of their stochastic traffic cellular au- 
tomaton (i.e., the STCA) pQ, led to a global adoptation 
of the TCA modeling scheme. One of the artefacts as- 
sociated with this STCA model, is that it gives rise to 
(many) unstable traffic jams A possible approach to 
achieve stable traffic jams, is to reduce the outflow from 
such jams, which can be accomplished by implementing 
so-called slow-to- start behavior. The name is derived from 
the fact that vehicles exiting jam fronts are obliged to wait 
a small amount of time. One implementation is based on 
making the noise parameter dependent on the velocity of a 
vehicle, leading to the development of the velocity depen- 
dent randomization (VDR) TCA. This TCA model more- 
over exhibits metastability and hysteresis phenomena 3 . 

In this context, our paper addresses the VDR-TCA 
slow-to-start model whose parameters take on a specific 
set of extreme values. Even though this bears little di- 



rect relevance for the understanding of traffic flows, it will 
lead to an induced 'anomalous' behavior and complex sys- 
tem dynamics, resulting in four distinct emergent phases. 
We study these on the basis of the system's relation be- 
tween density and flow (i.e., the fundamental diagram) 
which exhibits a non-concave region where, in contrast to 
the properties of congested traffic flow, density waves are 
propagated forwards. We finally investigate the system's 
tempo-spatial evolution in each of these four phases. 

Although some related studies exist (e.g., |415l6l7Kfl 
E]), the special behavior discussed in this paper has not 
been reported as such, as most previously done research 
is largely devoted to empirical and analytical discussions 
about these TCA models, operating under normal condi- 
tions. To strengthen our claims, we compare our results 
with the existing literature at the end of this paper. 

In section |5J we briefly describe the concept of a traf- 
fic cellular automaton and the experimental setup used 
when performing the simulations. The selected VDR-TCA 
model is presented in section [21 as well as an account 
of some behavioral characteristics under normal operat- 
ing conditions. The results of our various experiments 
with more complex system dynamics are subsequently pre- 
sented and extensively discussed in section 0] after which 
the paper concludes with a comparison with existing lit- 
erature in section and a summary in section EJ 
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2 Experimental setup 

In this section, we introduce the operational character- 
istics of a standard single-lane traffic cellular automaton 
model. To avoid confusion with some of the notations in 
existing literature, we explicitly state our definitions. 

2.1 Geometrical description 

Let us describe the operation of a single-lane traffic cellu- 
lar automaton as depicted in Fig.^ We assume N vehicles 
are driving on a circular lattice containing K cells, i.e., pe- 
riodic boundary conditions (each cell can be occupied by 
at most one vehicle at a time). Time and space are dis- 
cretized, with AT — 1 s and AX = 7.5 m, leading to a 
velocity discretization of AV = 27 km/h. Furthermore, 
the velocity Vi of a vehicle i is constrained to an integer 
in the range {0, . . . , Wmax}, with -Umax typically 5 cells/s 
(corresponding to 135 km/h). 
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Fig. 1. Schematic diagram of the operation of a single-lane 
traffic cellular automaton (TCA); here, the time axis is oriented 
downwards, the space axis extends to the right. The TCA's 
configuration is shown for two consecutive time steps t and 
t + during which two vehicles i and j propagate through the 
lattice. Without loss of generality, we denote the number of 
empty cells in front of vehicle i as its space gap g Si . 

Each vehicle i has a space headway h Si and a time 
headway h ti , defined as follows: 

h Si = Li + g Si , (1) 
hi = Pi+ 9u- (2) 

In these definitions, g Si and g^ denote the space and 
time gaps respectively; L t is the length of a vehicle and 
pi is the occupancy time of the vehicle (i.e., the time it 
'spends' in one cell). Note that in a traffic cellular automa- 
ton the space headway of a vehicle is always an integer 
number, representing a multiple of the spatial discretisa- 
tion AX in real world measurement units. So in a jam, it 
is taken to be equal to the space the vehicle occupies, i.e., 
h Si = 1 cell. 

Local interactions between individual vehicles in a traf- 
fic stream are modeled by means of a rule set. In this 
paper, we assume that all vehicles have the same physi- 
cal characteristics. The system's state is changed through 
synchronous position updates of all the vehicles, based on 
a rule set that reflects the car-following behavior. 

Most rule sets of TCA models do not use the space 
headway h Si or the space gap g Si , but are instead based 



on the number of empty cells di in front of a vehicle i. 
Keeping equation Q in mind, we therefore adopt the con- 
vention that, for a vehicle i its length Li = 1 cell. This 
means that when the vehicle is residing in a compact jam, 
its space headway h Si = 1 cell and its space gap is conse- 
quently g Si — cells. This abstraction gives us a rigorous 
justification to formulate the TCA's update rules more 
intuitively using space gaps. 

2.2 Performing measurements 

In order to characterize the behavior of a TCA model, 
we perform global measurements on the system's lattice. 
These measurements are expressed as macroscopic quan- 
tities, defining the global density k, the space mean speed 
v s , and the flow q as: 



i=l 

q = kv s . (5) 

The above measurements are calculated every time step, 
and they should be averaged over a large measurement 
period T gnll in order to allow the system to settle into an 
equilibrium. 

Correlation plots of these aggregate quantities lead 
to time-independent graphs conventionally called 'funda- 
mental diagrams'. And although we should more correctly 
refer to our measurements as points in a certain phase 
space (e.g., the (k,q) phase space), we will still use the 
terminology of 'fundamental diagram' in the remainder 
of this paper when we are in fact referring to this phase 
space. 

The previous global macroscopic measurements (den- 
sity, average speed, and flow) from equations © , Q , and 
(JSJ, can be related to the microscopic equations and 
(J2J) as follows: 

k oc h s , (6) 
q<xh t , (7) 

with h s and h t the average space and time headway re- 
spectively. Note that with respect to the time gaps and 
time headways, we will work in the remainder of this paper 
with the median instead of the arithmetic mean because 
the former gives more robust results when /i ti , g ti — > +oo 
for a vehicle i. 

All the fundamental diagrams in this paper, were cal- 
culated using systems of 10 3 cells. The first 10 3 s of each 
simulation were discarded in order to let initial transients 
die out; the system was then updated for T s j m = 10 4 s. 
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For a deeper insight into the behavior of the space 
mean speed v s , the average space gap ~g s , and the me- 
dian time gap ~g t , detailed histograms showing their dis- 
tributions are provided. These are interesting because in 
the existing literature (e.g., [TUIll 11112) 1 these distributions 
are only considered at several distinct global densities, 
whereas we show them for all densities. Each of our his- 
tograms is constructed by varying the global density k 
between 0.0 and 1.0, calculating the average speed, the 
average space gap and the median time gap for each sim- 
ulation run. A simulation run consists of 5 x 10 s (with a 
transient period of 500 s) on systems of 300 cells, varying 
the density in 150 steps. 

All the experiments were carried out with our Java 
software "Traffic Cellular Automata" 1 which can be found 
at http://smtca.dyns.cx |13|. 

3 Velocity dependent randomization 

In this section, the rule set of the VDR-TCA model is ex- 
plained, followed by a an overview of the model's tempo- 
spatial behavior and its related macroscopic quantities 
(i.e., the fundamental diagrams and the distributions of 
the speeds and the space and time gaps) under normal 
conditions. 



3.1 The VDR-TCA's rule set 

As indicated before, we focus our research on the VDR- 
TCA model for the implementation of the car-following 
behavior. The following equations (based on [3]) form its 
rule set; the rules are applied consecutively to all vehicles 
in parallel (i.e., synchronous updates): 

Rl: determine stochastic noise 

(vi(t- 1) = =^p' <-p , , s 

\vi(t-l)>Q=>>i/<-p, W 

R2: acceleration and braking 

Vi(t) <- mm{vi{t - 1) + l,g Si {t - l),v max }, (9) 

R3: randomization 

< P => «i(t) <- max{0,«i(t) - 1}, (10) 

R4: vehicle movement 

Xi{t) *-Xi(t-\)+13i{t). (11) 

In the above equations, Vi(t) is the speed of vehicle I 
at time t (i.e., in the current updated configuration of the 
system), Wmax is the maximum allowed speed, g Si denotes 
the space gap of vehicle i and Xi G {1, . . . , K} an integer 
number denoting its position in the lattice. In the third 
rule, equation (|TU)l . £At) G [0, 1[ denotes a uniform random 
number (specifically drawn for vehicle i at time t) and p' is 



the stochastic noise parameter, dependent on the vehicle's 
speed (po is called the slow-to-start probability and p the 
slowdown probability, with po,p G [0, 1]). 

In a nutshell, rule Rl, equation ||HJ|, determines the 
correct velocity dependent randomization. Rule R2, equa- 
tion © , states that a vehicle tries to increase its speed at 
each time step, as long as it hasn't reached its maximal 
speed and it has enough space headway. It also states that 
when a vehicle hasn't enough space headway, it abruptly 
adapts its speed in order to prevent a collision with the 
leading vehicle. The randomization parameter determined 
in equation JHJ, is now used in rule i?3, equation 1)10(1 . to 
introduce a stochastic component in the system: a vehicle 
will randomly slow down with probability p' . The last rule 
R4, equation 1)11)1 . isn't actually a 'real' rule; it just allows 
the vehicles to advance in the system. 

3.2 Normal behavioral characteristics 

Depending on their speed, vehicles are subject to different 
randomizations: typical metastable behavior results when 
Po 3> Pi meaning that stopped vehicles have to wait longer 
before they can continue their journey (i.e, they are 'slow- 
to-start'). This has the effect of a reduced outflow from a 
jam, so that, in a closed system, this leads to an equilib- 
rium and the formation of a compact jam. 

A 'capacity drop' takes place at the critical density, 
where traffic in its vicinity behaves in a metastable man- 
ner. This metastability is characterized by the fact that a 
sufficiently large disturbance of the fragile equilibrium can 
cause the flow to undergo a sudden decrease, correspond- 
ing to a first-order phase transition. The state of very high 
flow is then destroyed and the system settles into a phase 
separated state with a large megajam and a free-flow zone 
[3*lll4| . The large jam will persist as long as the density is 
not significantly lowered, meaning that recovery of traffic 
from congestion thus shows a hysteresis phenomenon |15j . 

If we look at the distribution of the vehicles' speeds, 
we get the histogram in Fig. Here we can clearly see 
the distinction between the free-flowing and the congested 
phase: the space mean speed remains constant at a high 
value, then encounters a sharp transition (i.e., the capac- 
ity drop), resulting in a steady declination as the global 
density increases. Once the compact jam is formed, the 
dominating speed quickly becomes zero (because vehicles 
are standing still inside the jam). 

Considering the distribution of the vehicles' space gaps, 
we get Fig. |31 because of their tight coupling in the VDR- 
TCA's rule set, the courses of both the space mean speed 
and the average space gaps are similar. Although the space 
gaps are rather large for low densities, at the critical den- 
sity k c they leave a small cluster around an optimal value 
of five cells. This corresponds to the necessary minimal 
space gap in order to travel at the maximum speed, avoid- 
ing a collision with the leader. An important observation 
is that no recorded space gaps exist between this cluster 
of five cells and the space gaps of zero cells inside the 
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Fig. 2. The distribution of the vehicles' speeds v, as a function of the global density k in the VDR-TCA (with po = 0.5 and 
p = 0.01). In the contourplot to the left, the thick solid line denotes the space mean speed, whereas the thin solid line shows its 
standard deviation. The grey regions denote the probability densities. The histograms (A) and (B) to the right, show two cross 
sections made in the left contourplot at k — 0.1325 and k = 0.4000 respectively: for example, in (A), the high concentration of 
probability mass at the histogram class v = 5 cells/s corresponds to the dark region in the upper left corner of the contourplot. 



compact jam. This means that there is a distinct phase 
separation taking place once beyond the critical density: 
vehicles are either completely in the free-flowing regime, 
or they are in the compact jam. 
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Fig. 3. The distribution of the vehicles' space gaps g s , as a 
function of the global density k in the VDR-TCA (with po = 
0.5 and p = 0.01). The thick solid line denotes the average of 
all the vehicles' space gaps, whereas the thin solid line shows 
its standard deviation. The grey regions denote the probability 
densities. 



Figure 0] shows the distribution of the vehicles' time 
gaps: when traffic is in the free-flowing regime, time gaps 
are high, but finite. As the critical density is approached, 
the median time gap first decreases (because the vehicles' 



speeds remain the same but their space gaps decrease). 
Once beyond the critical density, it increases towards in- 
finity because vehicles come to a full stop inside the com- 
pact jam. Just as with the space gaps, we can also ob- 
serve a small cluster around an optimal value. This op- 
timal time gap, is the time needed to travel the distance 
formed by the optimal space gap, at the maximum speed. 
This means that ~g t = g s -j-i?max = 1 s. Because the VDR- 
TCA incorporates stochastic noise, the median optimal 
time gap lies somewhat above the previously calculated 
value (g t « 1.2 s). 
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Fig. 4. The distribution of the vehicles' time gaps gt, as a 
function of the global density k in the VDR-TCA (with po — 
0.5 and p — 0.01). The thick solid line denotes the median of all 
the vehicles' time gaps. The grey regions denote the probability 
densities. 
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4 More complex system dynamics 

Most of the previous research dealt with the study of the 
behavioral characteristics of the VDR-TCA model oper- 
ating under normal conditions. We now turn our attention 
to the specific case in which the model's parameters take 
on extreme values po <C p, more specifically considering 
the limiting case where po = 0.0 and p — 1.0. 

We will first look at the change in tempo-spatial be- 
havior when p is increased towards 1.0, at which point a 
peculiar behavior is established in the system. We study 
the qualitative effects that the VDR-TCA's rule set has 
on individual vehicles, and discuss shortly the prevailing 
initial conditions. This is followed by a quantitative analy- 
sis using the (fc,g) fundamental diagram, leading us to the 
discovery of four distinct phases, having a non-concave re- 
gion with forward propagating density waves. More elabo- 
rate explanations are given based on the histograms of the 
space mean speed, the average space gap, and the median 
time gap. The section concludes with an analysis of the 
observations of the tempo-spatial behavior of the system 
in each of the four different phases (i.e., traffic regimes). 



4.1 Increasing the stochastic noise p 

Let us first consider the case in which po = 0-0 and where 
we vary p between 0.0 and 1.0. Figure [5] shows a time- 
space diagram where we simulated a system consisting of 
300 cells. As time advances (over a period of 580 s), the 
slowdown probability p is steadily increased from 0.0 to 1.0 
(the global density k was set to 0.1667 which is slightly be- 
low the critical density for a system with stochastic noise) . 



Time (580 s) with p G [0 -> 1] 

Fig. 5. A time-space diagram showing a system with a lattice 
of 300 cells (corresponding to 2.25 km); the visible time horizon 
is 580 seconds. As time advances, the slowdown probability p 
varies between 0.0 and 1.0 (the slow-to-start probability po was 
fixed at 0.0). The system's global density k is 0.1667. 

When p is very low, vehicles can keep driving in the 
free-flowing regime. As p is increased, small unstable jams 
occur. Increasing p even further, leads to even more pro- 
nounced jams. An important observation is that the prop- 
agation speed of these backward moving congestion waves 



increases; when p reaches approximately 0.5, their speed 
equals 0, meaning that jams stay fixed at a position. Note 
that as these jams are unstable, they can be created or dis- 
solved at arbitrary locations. Finally, as p tends to 1.0, we 
can see the emergence of forward propagating congestion 
waves; 'density' is now being carried in the direction of the 
traffic flow. These forward propagating waves form 'mov- 
ing blockades' that trap vehicles, which in turn leads to 
tightly packed clusters of vehicles that move steadily (but 
at a much slower pace than in the free-flowing regime). 
Note that this 'clustering' behavior is different from pla- 
tooning, which typically occurs when vehicles are driving 
close to each other at relatively high speeds 16 . In our 
case, the clusters of vehicles advance more slowly. 

4.2 Qualitative effects of the rule set 

From now on, we only consider the limiting case where 
Po = 0.0 andp =1.0. Let us now study the influence of the 
VDR-TCA's rule set on an individual vehicle i. Assuming 
Vi G {0, . . . , Umax}, the rules described in section l3~Tl lead 
to the following four general cases: 

- Case (1) with Vi{t - 1) = and g Si (t - 1) = 

In this case, the vehicle is residing inside a jam and 
rule R2 (acceleration and braking) plays a dominant 
part: the vehicle's speed Vi(t) remains 0. 

- Case (2) with v l (t - 1) = and g Sz (i - 1) > 

This situation may arise when, for example, a vehicle 
is at a jam's front. According to rule Rl, the stochastic 
noise parameter p' becomes 0.0. Rule R2 then results 
in an updated speed Vi(t) = min{l, g Si (t— 1)} = 1. Be- 
cause there is no stochastic noise present in this case, 
rule i?3 does not apply and the vehicle always advances 
one cell. 

- Case (3) with Vi(t- 1) > and g Sz (t - 1) = 
Because there is no space in front of the vehicle, it has 
to brake in order to avoid a collision. Rule R2 conse- 
quently abruptly decreases the vehicle's speed Vi(t) to 
0, as the vehicle needs to stop. 

- Case (4) with Vi(t - 1) > and g Sz (t - 1) > 

This case deserves special attention, as there are now 
two discriminating possibilities: 

- Case (4a) with v t (t - 1) < g Si (t - 1) 

According to rule Rl, the stochastic noise p' be- 
comes 1.0. Because the vehicle's speed is strictly 
less than its space gap, rule R2 becomes Vi(t) <— 
mm{vi(t — 1) + 1, Pmax}- Finally, rule R3 is applied 
which always reduces the speed calculated in rule 
R2 (constrained to 0). In order to understand what 
is happening, consider the speeds Vi(t — 1) and Vi(t) 
in the following table: 
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Vi(t - 



Vi(t) 



"max 
"max 
"max 

2 
1 



"max 
"max 
"max 

2 
1 



4.4.1 Effects on the fundamental diagrams 



We can clearly see that the maximum speed a vehi- 
cle can travel at, is constrained by Umax - 1, which 
corresponds to Umax — p' E] ■ From the table it 
follows that all vehicles traveling at Vi < Umax can 
neither accelerate nor decelerate: the vehicles' cur- 
rent speed is kept. 

- Case (4b) with v,(t - 1) > g Si {t - 1) 

Just as in the previous case (4a), the stochastic 
noise p' becomes 1.0. Rule R2 now changes to Vi(t) *- 
g Si (t — 1). Because rule R3 is always applied, this 
results in Vi(t) actually becoming g Si (t — 1) — 1 in- 
stead of just g Si (t — 1). So a vehicle always slows 
down too much (as opposed to solely avoiding a 
collision with its leader). 

Conclusion: considering the previously discussed four 
general cases, the most striking feature is that, accord- 
ing to case (4a), in a VDR-TCA model with po — 0.0 and 
p = 1.0, a moving vehicle can never increase its speed, 
i.e., v l (t) < Vi(t - 1). 



4.3 Effects of the initial conditions 

As already mentioned, in this paper we study the limiting 
case where p = 0.0 and p = 1.0. This case is special, in 
the sense that the system's behavior for light densities is 
extremely dependent on the initial conditions. 

After the system has settled into an equilibrium, the 
resulting flow is limited by the fact that the maximum 
speed of a vehicle in the system is always equal to Umax - 1< 
As proven in section 14.21 vehicles can never accelerate, 
which means that the slowest car in the system determines 
the maximum possible flow. Therefore, if the system (with 
Umax = 5 cells/s) is initialized with a homogeneous dis- 
tribution of the vehicles, then all of them will travel at 
Umax — 1 = 4 cells/s. If the system is initialized randomly, 
all vehicles will now travel at vL-.. 



1 cells/s with v' n 



the speed of the slowest vehicle in the system. 



4.4 Quantitative analysis 

Whereas the previous sections dealt with the effects of the 
VDR-TCA's changed rule set and the role of the initial 
conditions, this section considers the effects on the (k,q) 
fundamental diagram (i.e., three phase transitions and a 
non-concave region), as well as on the histograms of the 
space mean speed, the average space gap and the median 
time gap. 



Most existing (k,q) fundamental diagrams related to traf- 
fic flow models, show a concave course (although some 
(pedagogic) counter-examples such as |4I5I6I7I8I5| exist). 
Non-concavity (i.e., convexity) of a function / is defined 
as Vxi,X2 G dom/ | + x 2 )) < \f((xi) + f(x 2 ))- 

The property of concavity also holds true for the VDR- 
TCA model operating under normal conditions. However, 
when po <C p, this is no longer the case: in the limit where 
Po — 0.0 and p = 1.0, the fundamental diagram exhibits 
two distinct sharp peaks. Between these peaks, there ex- 
ists a region (77) + (777) where the fundamental diagram 
has a convex shape, as can be seen from Fig. [S] In region 
(III), traffic flow has a tendency to improve with increas- 
ing density. Note that we ignore the hysteretic behavior 
of the fundamental diagrams, as it gives no relevant con- 
tribution to the results of our experiments. 
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Fig. 6. The (k,q) fundamental diagrams of the VDR-TCA 
with a fixed slow-to-start probability po = 0.0 (the slowdown 
probability p is increased from 0.0 to 1.0). In the limiting case 
where p = 1.0, four distinct density regions (I)— (IV) appear. 



As the slowdown probability p is increased from 0.0 
to 1.0, the critical density - at which the transition from 
the free-flowing regime occurs - is shifted to lower val- 
ues (note that the magnitude of the capacity drop also 
diminishes). For a global density k — i, we can see that 
the speed of the backward propagating congestion waves 
increases, just as was visible in the time-space diagram 
of Fig. The speed of these characteristics is defined as 
dq/dk (i.e., the tangent to the fundamental diagram), and 
when p « 0.5, the sign of the speed is reversed, leading to 
the earlier mentioned forward propagating density waves. 
Furthermore, as is apparent in Fig. [ft] for high densities 
there exists a region (IV) in which the flow q is only de- 
pendent on k and not on the stochastic noise p: all the 
measurements coincide on this heavily congested branch. 



Sven Maerivoet, Bart De Moor: Non-concave fundamental diagrams and phase transitions in the VDR-TCA 



7 



Increasing the slowdown probability p has also an effect 
on the space mean speed in the free-flowing regime. As 
already stated in section l^l this speed is equal to Umax - 
p |17| . The average maximum speed in the free- flowing 
regime shifts downwards, reaching a value of 4 cells/s when 
p = 1.0. From then on, there are two regions (7) and (777) 
where the average speed remains constant. Note that, to 
be precise, region (7) actually contains a small capacity 
drop at a very low density, but we ignore this effect, thus 
treating region (7) in an overall manner. 

4.4.2 Effects on the histograms 

The distribution of the vehicles' speeds is shown in Fig.0 
where we can clearly observe two 'probability plateaus'. 
In the first region (7), vehicles' speeds are highly concen- 
trated in a small region around v s = 1 cell/s. As the global 
density increases in region (77), the space mean speed 
declines until it reaches region (777) where the second 
plateau is met at v s — 0.5 cell/s. From then on, it steadily 
decreases, reaching zero at the jam density (kj = 1.0). 
Note that in region (7), the standard deviation is zero, 
whereas it is non-zero but constant in region (III). This 
means that vehicles in the former region all drive at the 
same speed; in the latter region they drive at speeds alter- 
nating between and 1 cell/s (i.e., stop-and-go traffic). 



vehicles slowing down, because vehicles can never acceler- 
ate. 

Under normal operating conditions (i.e., po ^> p), a 
vehicle's average speed and average space gap show a 
high correlation in the congested density region beyond 
the critical density. This can be seen from the similar- 
ity between the histogram curves in Figs. [21 and [21 In high 
contrast with this, is the distribution of the vehicles' space 
gaps as in Fig. [SJ which shows a different scenario. More 
or less similar to the vehicles' speeds, we can observe the 
formation of three plateaus of constant space gaps in cer- 
tain density regions. These plateaus are located at 2, 1 and 
cells for density regions (7), (777), and (IV) respectively. 
Note that the average space gaps themselves are not con- 
stant in these regions, as opposed to the space mean speed. 
This is due to the fact that vehicles encounter waves of 
stop-and-go traffic, whereby the frequency of these waves 
increases as the global density is augmented. 

Another observation that we can make, is that the 
standard deviation goes to zero at the transition point 
between density regions (7) and (77). This means that, 
as expected, the traffic flow at this point consists of com- 
pletely homogeneous traffic, in which all the vehicles drive 
with the same space gap g s = 2 cells (and as already men- 
tioned, with the same speed v = 1 cell/s). 



1 



£ 3- 



0) p- 



0.0 0.1 0.2 



0.3 0.4 0.5 0.6 0.7 
Global density k 



0.9 1.0 



Fig. 7. The distribution of the vehicles' speeds v, as a function 
of the global density k in the VDR-TCA (with p = 0.0 and 
p — 1.0). The thick solid line denotes the space mean speed, 
whereas the thin solid line shows its standard deviation. The 
grey regions denote the probability densities. 




0.3 0.4 0.5 0.6 0.7 

Global density k 

Fig. 8. The distribution of the vehicles' space gaps g s , as a 
function of the global density k in the VDR-TCA (with po — 
0.0 and p — 1.0). The thick solid line denotes the average of 
all the vehicles' space gaps, whereas the thin solid line shows 
its standard deviation. The grey regions denote the probability 
densities. 



Considering the steep descending curve at the begin- 
ning of region (7), we state that although vehicles can 
drive at v = t>max — 1 = 4 cell/s under free-flowing condi- 
tions, they nonetheless all slow down as soon as at least 
one vehicle has a too small space gap. In other words, 
when g Si < t>max, case (4b) from section |4~2*1 applies and 
the vehicle slows down. This leads to a chain reaction of 



Just as with the previous histograms, there also ap- 
pear to be plateaus of concentrated probability mass in 
the distribution of the vehicles' time gaps in Fig. [5] As 
opposed to the standard behavior in Fig. the concen- 
tration in the first region (7) is more elongated and more 
or less completely flat. This is expected because the ma- 
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jority of the space mean speeds and the average space gaps 
remain constant in this region. Furthermore, once the hrst 
phase transition occurs, the median time gap g t — > +00 
as the space mean speeds and average space gaps tend to 
zero. This is expressed as the existence of a non-neglibible 
cluster of probability mass at the top of the histogram in 
Fig. 03 the concentration is formed by vehicles that are 
encountering the earlier mentioned stop-and-go waves, re- 
sulting in the fact that their time gaps periodically tend 
towards infinity. 
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Fig. 9. The distribution of the vehicles' time gaps gt, as a 
function of the global density k in the VDR-TCA (with po 
0,u and p = 1.0). The thick solid line denotes the median of all 
the vehicles' time gaps. The grey regions denote the probability 
densities. 




Fig. 10. Time-space diagram of the VDR-TCA model with 
po = 0.0 and p = 1.0. The shown lattice contains 300 cells 
(corresponding to 2.25 km), the visible time horizon is 580 
seconds. The density k — 0.2, corresponding to observations in 
density region (J), labeled free-flowing traffic (FFT). 



for the average space gap, explaining its rather 'nervous' 
behavior in Fig. |H1 

As the global density increases, the transition point 
between regions (J) and (II) is reached. At this point, 
each vehicle i has a speed Vi = 1 cell/s, and a space gap 
g Si = 2 cells. Because Vi < g Si , case (4a) from section l4~2*l 
applies. Using equations (JIJ and (JBJ, we can calculate the 
corresponding density: 



(J)-(JJ) 



= K l = (L- 



- 9s 



(12) 



Although the maximum speed any vehicle can (and 
will) travel at is limited to 1 cell/s, we still call this state 
'free-flowing traffic' (FFT) because no congestion waves 
are present in the system and none of the vehicles has to 
stop (see the time-space diagram in Fig. I1U|I . 



4.5 Typical tempo-spatial behavior 

Studying the (k,q) fundamental diagrams (Fig. EJ) in the 
previous section, we saw the emergence of four distinct 
density regions (I)-(IV) as p — ► 1.0 (formed by the A- 
symbols). In this section, we discuss the tempo-spatial 
properties that are intrinsic to these regions, relating them 
to the previously discussed histograms. 

As the global density k is increased, the system under- 
goes three consecutive phase transitions (between these 
four traffic regimes). The next four paragraphs detail the 
effects that appear in each density region (i.e., traffic re- 
gime), as well as the transitions that occur between them. 

4.5.1 Region (/) - free-flowing traffic [FFT] 

According to the histogram in Fig. |7| we can see that the 
speed of all the vehicles is the same (namely 1 cell/s) for 
moderately low densities. Although the standard devia- 
tion of the space mean speed is zero, this is not the case 



4.5.2 Region (//) - dilutely congested traffic [DCT] 

If we increase the density to k = -| + (i.e., adding one 
vehicle at the transition point), a new traffic regime is 
entered. In this regime, each extra vehicle leads to a back- 
ward propagating mini-jam of 3 stop-and-go cycles (see 
Fig. 1 1 1 p . bringing us to the description of 'dilutely con- 
gested traffic ' (DCT). 

In order to calculate the second transition point, we 
again observe the histogram of Fig. |SJ A vehicle's space 
gap now alternates between 1 and 2 cells in density re- 
gion (77). And because its speed alternates between and 

1 cell/s respectively, the vehicles' motions are controlled 
by cases (2) and (4a) from section l4~2l This means that 
stopped vehicles accelerate again in the next time step, 
after which they have to stop again, and so indefinitely 
repeating this cycle of stop-and-go behavior. 

Consider now a pair of adjacent driving vehicles i and 
j; it then follows from equation (JIJ that h s% = 1 + 1 = 

2 cells and h Sj =1 + 2 = 3 cells. So each pair of vehicles 
'occupies' 5 cells in the lattice, or 2.5 cells on average per 
vehicle. This leads to the second transition point being 
located at: 
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Fig. 11. Time-space diagram of the VDR-TCA model with 
po = 0.0 and p = 1.0. The shown lattice contains 300 cells 
(corresponding to 2.25 km), the visible time horizon is 580 
seconds. The density k = § + -j^, corresponding to observations 
in density region (//), labeled dilutely congested traffic (DCT). 



fc(HH(/H) = K 1 = — = 0.4. (13) 

As the density is increased towards fc(//)^(/77), the 
space mean speed decreases non-linearly. At the transition 
point itself, v s = 0.5 cells/s and the system is now com- 
pletely dominated by backward propagating dilute jams. 



4.5.3 Region (///) - densely advancing traffic [DAT] 

Adding one more vehicle at the transition point between 
regions (//) and (III), leads to surprising behavior: a for- 
ward moving jam emerges, traveling at a speed of 0.5 cells/s 
(see Fig. I12|l . Another artefact is that the cycle of alter- 
nating space gaps of 1 and 2 cells is broken, with the 
introduction of a zero space gap at the location of this 
new 'jam'. When the density reaches the third transition 
point, the system is completely filled with these forward 
moving 'jams' of dense traffic, leading to the description of 
'densely advancing traffic' (DAT). Because the available 
space is more optimally used by the vehicles, an increase 
of the density thus has a (temporary) beneficial effect on 
the global flow measured in the system. This kind of be- 
havior of forward moving density structures can also be 
observed in some models of anticipatory driving |16j . 

At the transition point itself, all vehicles exhibit the 
same behavior, comparable to the behavior at the sec- 
ond transition point. Traffic is more dense, as can be seen 
in the distribution of the space gaps in Fig. |H1 all ve- 
hicles have alternatingly g s = and g s = 1 cell (with 
corresponding speeds of and 1 cell/s respectively). One 
pair of adjacent driving vehicles i and j thus 'occupies' 
h Si + h Sj = (1 + 0) + (1 + 1) = 3 cells, or 1.5 cells on 
average per vehicle. The third transition point thus corre- 
sponds to: 

1 2 

k{ui)^{iv) = K = — = -. (14) 
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Fig. 12. Time-space diagram of the VDR-TCA model with 
po = 0.0 and p = 1.0. The shown lattice contains 300 cells 
(corresponding to 2.25 km), the visible time horizon is 580 
seconds. The density k = 0.4 + i, corresponding to observa- 
tions in density region (III), labeled densely advancing traffic 
(DAT). 

4.5.4 Region (IV) - heavily congested traffic [HCT] 

Finally, as the system's global density is pushed towards 
the jam density, each extra vehicle introduces at any point 
in time a backward propagating jam, consisting of a block 
of five consecutively stopped vehicles (see Fig. I13|) . The 
pattern of stable stop-and-go traffic gets destroyed, as ve- 
hicles remain stopped inside jams for longer time periods 
(cfr. density region (IV) in the distribution of the time 
gaps in Fig. 0, leading to the description of this traffic 
regime as 'heavily congested traffic' (HCT). 




Fig. 13. Time-space diagrams of the VDR-TCA model with 
po — 0.0 and p = 1.0. The shown lattice contains 300 cells 
(corresponding to 2.25 km), the visible time horizon is 580 
seconds. The density k = | + , corresponding to observations 
in density region (IV), labeled heavily congested traffic (HCT). 



5 Comparison with existing literature 

A final word should be said about the behavior of the 
four phases discussed in this paper. Although the study is 
based on a well-known traffic flow model (i.e., the VDR- 
TCA model), it seems that this behavior is drastically dif- 
ferent from that observed in real-life traffic. The research 
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elucidated in this paper, is therefore relevant as it can be 
compared to phase transitions in other types of granular 
media, described by cellular automata. 

Relating our findings to similar considerations in liter- 
ature, we note the following aspects: Grabolus [HJ gives a 
thorough discussion of several variants of the STCA, in- 
cluding so-called 'fast-to-start' TCA models that give rise 
to forward propagating density waves; the study refers to 
the two-fluid theory nfS| as a possible explanation of this 
phenomenon. Gray and Griffeath [§] discuss the origin of 
non-concavity in the (k,q) fundamental diagram, where 
they propose an explanation in which the rear end of a 
jam is unstable, whereas its front is stable and growing. 
Leveque |7| incorporates a non-concave fundamental di- 
agram itself in a macroscopic model, thereby resembling 
night driving. Nishinari et al [B] consider the tempo-spatial 
organization of ants using an ant trail model based on 
pheromones. In their model, the average speed of the ants 
varies non-monotonically with their density, leading to an 
inflection point in the {k,q) fundamental diagram. Zhang 
|!5] questions the property of anisotropy in multi-lane traf- 
fic, leading to strikingly similar non-concave (k,q) funda- 
mental diagrams. 

Finally, Awazu PI] investigated the flow of various com- 
plex particles using a simple meta-model based on Wol- 
fram's CA-rule 184, but extended with rules that govern 
the change in speed of individual particles. As a result, 
he classified three types of fundamental diagrams, namely 
two phases (2P), three phases (3P), and four phases (4P) 
type systems. Of these types, the 4P-type complex gran- 
ular particle systems closely resemble the regimes discov- 
ered in this paper. Awazu discusses the types of regimes 
and the transitions between them, using the terminology 
of 'dilute slugs' for the slow moving jams in the DCT- 
phase (see section 14. 5.2(1 . which he calls the 'dilute jam- 
flow state'. Analogously, there are 'advancing slugs' in the 
'advancing jam- flow state' (related to our DAT-phase in 
section B.5.811 and the fourth phase is called the 'hard jam- 
flow state' (related to our HCT-phase in section PL 5. 4fl . 

As stated before, the behavior here is different from 
that of real-life traffic flows. Awazu expects that these 
transitions appear if there are more complex interactions 
between the particles. In the 4P-type systems, attractive 
forces between close neighbouring particles and an effec- 
tive resistance on moving neighbouring particles, lead to 
the realization of the previously discussed flow phases. 



6 Summary 

In this paper, we first showed the behavioral character- 
istics resulting from the velocity dependent randomiza- 
tion traffic cellular automaton (VDR-TCA) model, oper- 
ating under normal conditions (i.e., po 3> p). Then we in- 
vestigated the more complex system dynamics that arise 
from this traffic flow model in the exceptional cases when 
Po *C P- This behavior was quantitavely compared against 



the VDR-TCA's normal operation, using classical funda- 
mental diagrams and histograms that show the distribu- 
tion of the speeds, space, and time gaps. 

Our main investigations were primarily directed at the 
limiting case where po = 0.0 and p = 1.0. We discov- 
ered the emergence of four different traffic regimes. These 
regimes were individually studied using diagrams that show 
the evolution of their tempo-spatial behavioral character- 
istics. This resulted in the following classification: free- 
flowing traffic (FFT), dilutely congested traffic (DCT), 
densely advancing traffic (DAT), and heavily congested 
traffic (HCT). Our main conclusions here are: 

— all four phases share the common property that moving 
vehicles can never increase their speed once the system 
has settled into an equilibrium, 

— the DAT regime experiences forward propagating den- 
sity waves, corresponding to a non-concave region in 
the system's flow-density relation. 

Comparing our results with those in existing literature, 
we conclude that the work of Awazu Pjj , dealing with a cel- 
lular automaton model of various complex particles, gives 
the closest resemblance to our four phases. 

Note that a more detailed account of our research (in- 
cluding (k,v s ) fundamental diagrams, more on the influ- 
ence of the maximum speed and the use of a suitable order 
parameter to track the phase transitions) can be found in 
our technical report |19j . 
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